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Motivated by the generalization of quantum theory for the case of 
non-Herniitian Hamiltonians with PT symmetry, we show how a clas- 
^^ ■ sical cosmological model describes a smooth transition from ordinary 

^D I dark energy to the phantom one. The model is based on a classical 

^^ ■ complex Lagrangian of a scalar held. Specihc symmetry properties 

analogous to PT in non-Hermitian quantum mechanics lead to purely 



o 



O ■ real equation of motion. 
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c^ ■ 1 Introduction 



Complex (non-Hermitian) Hamiltonians with PT symmetry have been vig- 
orously investigated in quantum mechanics and quantum field theory pP . A 
possibility of applications to quantum cosmology have been pointed out in 
P] . In the present contribution we focus attention on complex classical field 
theory. We explore the use of a particular complex scalar field Lagrangian, 
which has real solutions of the classical equations of motion. Thereby we 
provide a cosmological model describing in a natural way an evolution from 
the Big Bang to the Big Rip involving the transition from normal matter to 
phantom matter, crossing smoothly the phantom divide line. 

The interest of our approach is related to its focusing on the intersection 
between two important fields of research. The basic idea can be introduced 



qualitatively as follows: given a complex Lagrangian of complex scalar field, 
with the complex potential V{(p), if V{(f)) is PT symmetric in the sense that 
V{x) is PT symmetric in quantum mechanics, e.g. V{(p) ~ exp{ia(p) with a 
real, then Lagrangian becomes real for purely imaginary (f> and furthermore 
the kinetic term acquires a negative sign. In order to avoid any confusion 
we stress that our arguments are qualitative and PT symmetry is not to 
be intended in a literal way since we will deal with spatially homogeneous 
scalar fields and we have introduced these notations to refer the reader to 
the type of potentials and to their symmetry properties we will consider. 

The framework in which we work is the general relativity and classical 
cosmology. We will provide a self-contained introduction to these topics in 
order that an average reader can follow the presentation without too many 
difficulties. 

2 Complex Lagrangians in classical field theory 
and cosmology 

Let us consider a non-Hermitian (complex) Lagrangian of a scalar field 
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with the corresponding action, 



L = -d^4>d>'4>*-v{<p,^*), 



S{cP,cP\g) = f d^x.f^\{L + ^R{g)) 
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where ||5f|| stands for the determinant of a metric g^'^ and R{g) is the scalar 
curvature term and the Newton gravitational constant is normalized to S/Svr 
to simplify the Friedmann equations further on. 

We employ potentials V{<l),(j)*) satisfying the invariance condition 

(y(</>,0*)r = y(0*,</.), (3) 

while the condition 

(y(0,0*)r = y(0,0*), (4) 

is not satisfied. For example, such potential can have a form 

v{(t>,^*) = Vi{4> + <P*)V2{4>-<p*), (5) 

where Vi and V2 are real functions of their arguments. If one defines 

<Pi^^{<P + 4>*), (6) 



and 

02^^(0 -</>*), (7) 

one can consider potentials of the form 

y((^,(/.*) = yo(0i)exp(ia02), (8) 

where q is a real parameter. In the last equation one can recognize the link 
to the so called PT symmetric potentials. 

Here, the functions (j)i and (f)2 appear as the real and the imaginary parts 
of the complex scalar field (j), however, in what follows, we shall treat them 
as independent spatially homogeneous complex variables depending only on 
the time parameter t. 

We shall consider a flat spatially homogeneous Friedmann universe with 
the metric 

ds^ = dt^ - a^{t)df, (9) 

satisfying the Friedmann equation 

/i' = ^'/'? + ^02 + ^o(</'i)exp(ia(/.2). (10) 

Here the variable a{t) represents a cosmological radius of the universe and 
the Hubble variable h{t) = ^ characterizes the velocity of expansion of the 
universe. The Friedmann equation (jlUj) is nothing but Einstein equation, 
for the universe filled by scalar fields. 

The equations of motion for fields (/)i and 02 have the form 

<pi + 3h(pi + V^{(pi) exp(ia(/.2) = 0, (11) 

z(/>2 + 3ih(j)2 — aVo{(j)i) ex.p{iacj)2) = 0. (12) 

Eqs. (|1U|) . (|11|) and (|12j) are obtained by variation of the action Q with 
the lagrangian (^ and the potential ^ with respect to the metric, and the 
scalar field variables (pi and (/>2. 

3 Cosmological solutions for accelerated universes 

Let us notice, that the system of equations (fTT|) .((T ^ . HlU() can have a solution 
where (piit) is real, while the (p2 is imaginary, or, in other words 

Mt) = -im, (13) 



where ^(t) is a real function. In terms of these two functions, our system of 
two real equations can be rewritten as 



\^l - \e + ^o(</>i) expK)0i + ^o(</'i) exp(aO = 0, (14) 



I + ^\j\4>l - \e + Vo{4>i) exp(aOe - «l^o('/'i) expK) = 0. (15) 

Now, substituting (/>2(t) from Eq. (^SJ into Eq. (|l()j) we have the follow- 
ing expression for the energy density 

e = h' = ]^^l- i^2 + y^^^^) exp«). (16) 

Thus, we have paved the way to convert an action with a complex potential 
into the action with real potential and hyperbolic structure of the kinetic 
term. Notice that in the flat Friedmann universe the energy density is always 
positive, being proportional to squared Hubble variable which should be real 
because of reality of geometry. 

Let us emphasize that the solution which we are looking for, namely, the 
imaginary solution for 02 is the solution which can have sense for the whole 
system of equations of Klein-Gordon type (|1H) . (|12|) and that of Friedmann 
(jlUj) . because it makes the energy density real and positive. One can trace 
some kind of analogy with the non-Hermiatian PT symmetric quantum the- 
ory with real spectra bounded from below. Our motivation is just to combine 
the idea of i) providing physical application for "PT" -symmetric potentials 
and of ii) building of proper framework for conversion of the elliptic struc- 
ture of the kinetic term for the scalar field to the hyperbolic one. In fact, 
it is crucial to start from a complex Lagrangian with the above mentioned 
symmetry in order that Eqs. (|11() and H12|) become real after the rotation 

h{t) = -im- 

We shall see in the section 4 that one can construct such solutions so that 
the originally complex Lagrangian becomes real on classical configurations 
while one of the scalar components obtains the ghost (negative) sign of 
kinetic energy. Thereby we recover a more conventional phantom matter 
starting from the complex matter with normal kinetic energy. 

Now, coming back to the Einstein equations we should remember that 
the pressure will be equal 

P=^0?-^^'-Vb(</'i)exp(aO. (17) 



The pressure in cosmology can be negative. Moreover, the so called dark 
energy, responsible for recently discovered phenomenon of cosmic accelera- 
tion P], is characterized by negative pressure such that w = p/e < —1/3 
[1]. The value t(; = — 1 is nothing but the cosmological constant, while the 
dark matter with w < —1 was dubbed as phantom energy j2|. The phantom 
models have some unusual properties: to realize them one often uses the 
phantom scalar field with the negative sign of kinetic term; in many models 
the presence of the phantom dark energy implies the existence of the future 
Big Rip cosmological singularity [O]- A cosmological evolution where dark 
energy undergoes the transition from w > —1 to w < —1 implies some par- 
ticular properties of the corresponding field-theoretical model and is called 
crossing of phantom divide line [7j. 

From now on we show that the model introduced above is suitable for 
the description of the phenomenon of the phantom divide line crossing (for 
further detail see [H]). 

4 Crossing of the phantom divide Hne 

It is easy to see that if 0f < ^^ the pressure will be negative and p/e < 
— 1, satisfying the phantom equation of state. Instead, when (j)^ > ^^, the 
ratio between the pressure and energy density exceeds —1 and, hence, the 
condition 

0? = f (18) 

corresponds exactly to the phantom divide line, which can be crossed dy- 
namically during the evolution of the field components (j)i{t) and ^(t). 

We provide now a simple realization of this idea by an exactly solvable 
cosmological model by implementing the technique for construction of po- 
tentials for a given cosmological evolution [Sj. It is convenient to start with 
a cosmological evolution as given by the following expression for the Hubble 
variable: 

The evolution begins at t = 0, which represents a standard initial Big Bang 
cosmological singularity, and comes to an end in the Big Rip type singularity 
at t = tpi,. The derivative of the Hubble variable 

. _ A(2t - tr,) 



vanishes at 

tp = Y (21) 

when the universe crosses the phantom divide hne. 

Next, we can write down the standard formulae connecting the energy 
density and the pressure to the Hubble variable and its time derivative: 

f4 + W0.-' = ^^ = J,iJ^. (22) 

2 2 "^^^ 3 3t^{tR-tf ^ ' 

The expression for the potential Vo(0i) follows 



The kinetic term satisfies the equation 



2yl(2t - tR ) 
'3t^itR-t) 



i' = - :?.:.% ■ (25) 



It is convenient to begin the construction with the solution for ^. Taking 
into account the formulae ()19() and H24|) Eq. ()15() can be rewritten as 

Introducing a new parameter 
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(27) 


Eq. pfi)) looks hke 
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= 0, 


(28) 
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where 




y^e 






(29) 



The solution of Eq. H28|) is given by 



_ amtRJtR - tr f ^^ (2t + (m - l)^^)^'"-^ 



Before considering the concrete values of m, notice that the equation of 
state parameter w in the vicinity of the initial Big Bang singularity behaves 
as 

w = -l + —, (31) 

m 

while approaching the final Big Rip singularity this parameter behaves as 

2 . ^ 

w = -l . 32 

m 

Notice that the range for w does not depend on a, depending only on the 
value of the parameter m, which relates the scales of the Hubble variable h 
and of the time of existence of the universe t^. 

Remarkably, an integral in the right-hand side of Eq. (|30jl is calculable 
analytically 

i = ^^ (33) 

while 

rvTD 

e = — (logt-log(tH-t)). (34) 

From now on the parameter t will be dimensionless. Inclusion of charac- 
teristic time does not change the structure of the potential because of its 
exponential dependence on ^. Substituting the expression (|55|) into Eq. (^^ 

one has 

•2 mtR{{a^m + l^)tR-mt) 

^^ = ,lt^itR-tf • ('') 

For the case a^m = 18 the function 0i(t) can be easily found from Eq. ()35|) 
and it looks like follows: 



^^ = ±l^/^ArctanhJ^^^. (36) 

One can choose the positive sign in Eq. H36() without loosing the generality. 
Inverting Eq. ()36j) we obtain the dependence of the time parameter as 
a function of 0i 

Ay/m 

Substituting expressions (|S7|) and (|51)) into Eq. (pi)) we can obtain the 
explicit expression for the potential Vb (</>!): 

2cosh6 ^ (2 + 17cosh2 ^) 
Fo(</)i) = ^^^^^-2 '-^^. (38) 



We would like to emphasize that this potential is real and even. It is inter- 
esting that the time dependence of (j)i{t) could be found also for an arbitrary 
value of the parameter m, but for a'^m > 18 this dependence cannot be re- 
versed analytically and, hence, one cannot obtain the explicit form of the 
potential Vo(0i). 

5 Concluding remarks and perspectives 

The phantom model building has involved many different ideas. Here we 
have presented a rather simple and natural cosmological toy model, linked to 
and inspired by such an intensively developing branch of quantum mechanics 
and quantum field theory as the study of non-Hermitian, but CPT (or PT) 
symmetric models. Notice that there is an analogy between the manner in 
which the complexity of the original Lagrangian with the standard kinetic 
term is transformed into phantom-like Langrangian, which is real but has 
a negative kinetic energy term and the equivalence between PT-symmetric 
quantum Hamiltonians and Hermitian Hamiltonian with variable effective 
mass (see e.g. [T0]1. 

In this paper we have focused essentially on the classical field theory with 
a complex potential satisfying the invariance property of Eq. (jSJ. Beyond 
the classical limit one could speculate on how quantum fluctuations 5(j) = 
r]i{x) + ir]2{x) may preserve the consistency of this approach. Assume that 
the fields of fluctuations respect the initial and final conditions on a classical 
solution, 7/1(0) = r?2(0) = r]i{tfi) = ri2{tfi) = . Then the second variation 
of the action Q reads, 

S^^\m,V2) = \j d''xa\t)Tfi-d^t-2,h{t)dt-V'^^\(t)i,cM2,c)\. (39) 

where r/-^ = (r/i(x),ry2(x)) is the transposed field , ct)i^c = 0i(i)) '^2,c = —i^{t) 
are classical solutions and the matrix V^"^' reads, 

The quadratic form (|39() is symmetric, with a non-Hermitian but pseudo- 
Hermitian matrix asV^'^'a^ = {V^'^')' and the latter fact makes it possible 
to get real eigenvalues (or pairs of complex conjugated ones) of the energy 
operator in (|39j) . In so far as this energy operator is symmetric and has a 



2x2 matrix form, one can diagonalize it^ with an (in general, non-local) 
orthogonal transformation O so that 00^ = I. 

Because the second variation of potential V^"^' has complex matrix ele- 
ments, the eigenvectors for a particular real eigenvalue will be also complex. 
The correct way to perform the variation is: 

first to make an appropriate complex deformation of the integration con- 
tour of variables r]i{x),rj2{x) — > fji{x),fJ2{x) so that the latter complex 
variables give rise to real ones r] = Ofj; 

next to perform the diagonalization and to end up with a well defined, 
real energy operator, hopefully with positive masses. 

To realize this program one has to solve the one-dimensional matrix 
Schrodinger-like equation. To give a concrete idea of what we have in mind 
we present a solvable example of potential (only mildly related to the previ- 
ous discussion), V = exp(\/l + a'^ (pi) exp(ia(/)2)- Its second variation matrix 
is constant up to an overall factor. 



1 + a^ iaVl + 



\ lay I + a^ -a^ / 

Its eigenvalues are 0, 1 and the normalized eigenvectors are 
e^ = {ia, — \/l + a"^), ej = (a/1 + a'^,ia) so that efcj = Sij . The diago- 
nalization is realized by the complex orthogonal matrix 



0=1 /^ -VTT^Y (42) 

\ vl + a la I 

The above mentioned contour is given by f/ = O rj for arbitrary real r]. Thus, 
in spite of non-Hermiticity of the matrix of second variation the relevant 
deformation of the contour makes the action real with positive kinetic terms 
and masses. 

As a last remark we would like to point out that in quantum mechanics, 
models with CPT symmetry have been recently introduced in paper [12^ with 
the difference that the charge operator is a differential operator contributing 
to a definition of a pseudometric operator, whereas in the present approach 
the "charge" conjugation is related to a kind of internal degree of freedom. 

Finally, let us recapitulate various steps of our approach. 
1. We relate the possibility of a (de)- phantomization to a rotation of one of 

^ While it is not true in general for complex symmetric matrix in arbitrary dimension 
as discussed for example in |11|. in the case of a quadratic form in Ea. H39|l . one can 
straightforwardly show the above statement by explicit construction of the corresponding 
complex orthogonal 2x2 matrix. 



the compenent of the scalar field, which gives to the kinetic term a hyperbolic 
structure. 

2. The preceding point provides compelling reasons to start from complex 
PT symmetric Lagrangian. 

3. We succeed to reduce the problem to two coupled scalar fields: of course 
one could have started directly from these coupled equations, but in this 
case, underlying CPT symmetry and CP breaking would be hidden 

We do not claim uniqueness of our approach (see e.g. the recent paper [T^ . 
based on the superalgebraic approach and using Grassmann vector fields as 
well as a more general complexification involving the space-time coordinates 
[llj). but a definite consequence turns out to be the relation between possible 
phantomization and CP-breaking. 
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